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Spectral features of solar wind turbulent plasma 
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ABSTRACT 

Spectral properties of a fully compressible solar wind Hall Magnetohydrodynamic 
plasma are investigated by means of time dependent three dimensional Hall MHD 
simulations. Our simulations, in agreement with spacecraft data, identify a spectral 
break in turbulence spectra at characteristic length-scales associated with electromag- 
netic fluctuations that are smaller than the ion gyroradius. In this regime, our 3D 
simulations show that turbulent spectral cascades in the presence of a mean magnetic 
field follow an omnidirectional anisotropic inertial range spectrum close to /c^^/'^. The 
onset of the spectral break in our simulations can be ascribed to the presence of nonlin- 
ear Hall interactions that modify the spectral cascades. Our simulations further show 
that the underlying charachteristic turbulent fluctuations are spectrally anisotropic, 
the extent of which depends critically on the local wavenumber. The fluctuations asso- 
ciated with length scales smaller than the ion gyroradius are highly compressible and 
tend to exhibit a near equipartition in the velocity and magnetic fields. Finally, we find 
that the orientation of velocity and magnetic field fluctuations critically determine the 
character of nonlinear interactions that predominantly govern a Hall MHD plasma, 
like the solar wind. 

Key words: (magnetohydrodynamics) MHD, (Sun:) solar wind. Sun: magnetic fields, 
ISM: magnetic fields 



1 INTRODUCTION 

The solar wind plasma is predominantly in a turbulent state. 
Nonlinear turbulent processes in the magnetized solar wind 
plasma plasma fluid yield a multitude of spatial and tempo- 
ral length-scales associated with an admixture of waves, fluc- 
tuations, structures and nonlinear turbulent interactions. In- 
situ spacecraft measurements (Matthaeus & Brown 1988, 
Goldstein et al 1995, Ghosh et al 1996) reveal that the so- 
lar wind fluctuations, extending over several orders of mag- 
nitude in frequency and wavenumber, when described by 
power spectral density (PSD) spectrum, can be divided into 
three distinct regions (Goldstein et al 1995, Leamon et al 
1999) depending on the frequency and wavenumber. The 
first region corresponds to a flatter spectrum, associated 
with lower frequencies, and it is consistent with (where 
k is wavenumber). Second region follows and extends to the 
ion/proton gyrofrequency, and the spectral slope has an in- 
dex ranging from -3/2 to -5/3. This region is typically char- 
acterized as corresponding to fully developed turbulence, 
and can be described by the usual incompressible magne- 
tohydrodynamic (MHD) description. The turbulent interac- 
tions in this regime are governed entirely by Alfvenic cas- 
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cades. Spacecraft observations (Leamon et al 1999, Bale et 
al. 2005, Alexandrova et al 2007, Sahraoui et al 2009) fur- 
ther reveal that at length scales beyond the MHD regime, 
i.e. length scales less than ion gyro radius [kpi ^ 1) and 
temporal scales greater than the ion cyclotron frequency 
Lo > Uci = eBo/nieC, (where k, pi,ujci,e, Bo,me,c are re- 
spectively characteristic mode, ion gyroradius, ion cyclotron 
frequency, electronic charge, mean magnetic field, mass of 
electron, and speed of light), the spectrum exhibits a spec- 
tral break, and the spectral index of the solar wind tur- 
bulent fluctuations varies between -2 and -5 (Smith et al 
1990, Goldstein et al 1994, Leamon et al 1999, Bale et al. 
2005, Shaikh & Shukla 2009, Sahraoui et al 2009). Higher 
time resolution observations flnd that at the spectral break, 
Alfvenic MHD cascades (Smith et al 1990, Goldstein et al 
1994, Leamon et al 1999) close. The characteristic modes in 
this region appear to evolve typically on timescales associ- 
ated with dispersive kinetic Alfvenic fluctuations. 

The onset of the second or the kinetic Alfven iner- 
tial range is not understood. Some suggestions have how- 
ever been made. The spectral break may result from en- 
ergy transfer processes associated with possibly kinetic 
Alfven waves (KAWs) (Hasegawa 1976), electromagnetic 
ion-cyclotron- Alfy en (EMICA) waves ()Garv et al 20081 : 
IWu fc Yoon 20071 ). or by fluctuations described by a Hall 
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MHD (HMHD) plasma model (Alexandrova et al 2007, 
2008; Shaikh & Shukla 2008, 2008a). Stawicki et al (2001) 
argue that Alfven fluctuations are suppressed by pro- 
ton cyclotron damping at intermediate wavenumbers so 
the observed power spectra are likely to comprise weakly 
damped dispersive magnetosonic and/or whistler waves (un- 
like Alfven waves). Beinroth & Neubauer (1981) and Den- 
skat & Neubauer (1982) have reported the presence of 
whistler waves based on Ifelios 1 & 2 observations in this 
high frequency regime. A comprehensive data analysis by 
Goldstein et al. (1994), based on correlations of sign of mag- 
netic helicity with direction of magnetic field, indicates the 
possibility of the existence of multiscale waves (Alfvenic, 
whistlers and cyclotron waves) with a single polarization in 
the dissipation regime. Counter intuitively, in the u < uici 
regime, or Alfvenic regime, Ifowes et al. (2008) noted the 
possibility that highly obliquely propagating KAWs are 
present (with uj ^ ujd) making the possibility that damping 
of ion cyclotron waves is responsible for the spectral break- 
point questionable. 

Fluid (Shaikh & Shukla 2009) and kinetic (Ifowes et 
al. 2008) simulations, in qualitative agreement with space- 
craft data described as above, have been able to obtain 
the spectral break point near the characteristic turbulent 
length scales that are comparable with ion inertial length 
scale (di). These simulations demonstrated a Kolmogorov- 
like k~^^^ spectra for the length scales larger than the ion 
inertial length scales where MHD is typically a valid de- 
scription. By contrast, smaller (than di) scales are shown 
to follow a steeper spectrum that is close to k~''^^ (Howes 
et al. 2008, Shaikh & Shukla 2009). Spacecraft data and 
simulations thus reveal that migration of turbulent energy 
proceeds essentially through different regions in the fc-space, 
i.e. k~^,k~^^^ and k~'^^^. Needless to mention that turbu- 
lent cascade does not entirely terminate immediately beyond 
the k~'^^^ spectrum. Fluid and kinetic simulations (Biskamp 
1996, Gaftier 2006, Gaftier & Buchlin 2007, Cho and Lazar- 
ian 2004, Shaikh & Zank 2005, Shaikh 2009, Gary et al. 2008, 
Saito et al. 2008, Howes et al. 2008) have further shown that 
spectral transfer of energy extends even beyond the k'"^^^ 
spectrum and it is governed predominantly by small scale, 
high frequency, whistler turbulence. The latter also exhibits 
a definite power law. 

The physical processes describing MHD, KAW or Hall 
MHD and whistler spectra are rich and complex. They differ 
significantly from others and continue to pose serious chal- 
lenges in our understanding of multiscale solar wind turbu- 
lence. One of the major goals of this paper is to describe 
the connection between difi'erent scales associated with the 
MHD, KAW or Hall MHD and whistler spectra. In the fol- 
lowing, we describe extended part of the spectra that are 
predicted by theory and simulations. 



2 EXTENDED COMPOSITE SPECTRA 

Theory and simulations indicate that turbulent fluctuations 
in the high frequency and kpi ^ 1 (where pi is ion gyro ra- 
dius) regime correspond to a regime in which electron mo- 
tions are decoupled from the ion motions (Kingsep et al 
1990, Biskamp et al 1996, Dastgeer et al 2000a, Dastgeer et 
al 2000b, Shaikh & Zank 2003, Cho & Lazarian 2004, Saito 



Figure 1. Schematic of the power spectral density (PSD) com- 
posite spectrum in the solar wind turbulent plasma as a function 
of frequency (wavenumber). Several distinct regions are identi- 
fied with what are thought to be the dominant energy transfer 
mechanism for that particular region. The nonlinear processes 
associated with the transition from region II (MHD regime) to 
region III (kinetic or Hall MHD regime) are not yet fully under- 
stood. The power spectra in region III vary from to 
The boundary of region III and IV identifies where electron and 
ion motions are decoupled. Regions IV and V are identified as 
whistler cascade regimes. 

et al 2008, Gary et al 2008). Correspondingly, ions are es- 
sentially unmagnetized and can be treated as an immobile 
neutralizing background fluid. This regime corresponds to 
the whistler wave band of the spectrum and comprises char- 
acteristic scales that are smaller than those that describe 
MHD, KAW or Hall MHD processes. An extended compos- 
ite schematic thus describing the whistler modes spectra, in 
addition to the observed k~^,k~^^^ and k~''^^ spectra, is 
shown in Fig. ((T}. Speciflcally, regions IV and V in Fig. ((TJ 
identify characteristic modes that are relevant for the de- 
scription of whistler wave turbulence (Biskamp et al 1996, 
Shaikh & Zank 2005, Shaikh 2009). By contrast, regions I, 
II and HI describe respectively k~^,k~^^^ and k~'^^^ spec- 
tra and are consistent with the observations (Leamon et al 
1999, Bale et al. 2005, Alexandrova et al 2007, Sahraoui et 
al 2009). The boundary of regions HI and IV represents a 
wavenumber band in spectral space that corresponds to the 
decoupling of electron and ion motions. The wavenumbers 
above this boundary characterize the onset of whistler tur- 
bulence. The spectral cascades associated with whistler tur- 
bulence are described extensively by Biskamp et al (1996), 
Dastgeer et al (2000a), Dastgeer et al (2000b), Shaikh & 
Zank (2003), Shaikh & Zank (2005), Shaikh (2009a), Shaikh 
(20009b). Cho & Lazarian 2004 describe scale dependent 
anisotropy that is mediated by whistler waves in the con- 
text of electron MHD plasma. Gary et al. (2008) and Saito 
et al. (2008) have reported two-dimensional electromag- 
netic particle-in-cell simulations of electron MHD model 
to demonstrate the forward cascade of whistler turbulence. 
Their work show that magnetic spectra of the cascading 
fluctuations become more anisotropic with increasing fine- 
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tuation energy. Interestingly, whistler turbulence associated 
with longer wavelengths in region IV exhibits a power spec- 
trum k~'^^^ that is similar to the short wavelength spectrum 
of kinetic Alfven waves (KAW), as shown in region III of Fig. 
([l]). The underlying physical processes, responsible for the 
spectrum differ significantly for KAW and whistler waves. 
These differences are discussed in Section 7. 

It is to be noted that the Hall MHD description of mag- 
netized plasma is valid upto region III where characteristic 
turbulent scales are smaller than ion inertial length scales 
{kdi > 1). Beyond this point, the high frequency motion of 
plasma is governed predominantly by the electron motions 
only. By contrast, ions form static neutralizing background. 
Consequently, the ion motions decouple significantly from 
the electrons. These aspects of the spectra, depicted neces- 
sarily by regions IV & V in Fig (1), can be described ade- 
quately by whistler wave model. The Hall MHD models are 
therefore not applicable in regions IV, V and beyond. Nei- 
ther they can describe kinetic physics associated with the 
dissipative regime. Since the high frequency regime (i.e. re- 
gions IV & V) is dominated by the electron motions, there 
exists intrinsic length scale corresponding to electron iner- 
ital length scale = c/cup^ (where c is the speed of light 
and u>p^ is the electron plasma frequency). The characteris- 
tic turbulent length scales in regions IV & V are essentially 
comparable with and therefore they can describe scales 
larger (i.e. kd^ < 1 in region IV) and smaller (i.e. kde > 1 
in region V) than the electron inertial scale. While whistler 
wave model can describe nonlinear processes associated with 
length scales as small as the electron inertial length scale, 
they fail to describe finite electron Larmour radii effects for 
which a fully kinetic description of plasma must be sought. 

The schematic of solar wind turbulence depicted in Fig. 
lU raises numerous unresolved questions. Beside those de- 
scribed above, we do not understand what leads to the de- 
coupling of ion and electron motions near the boundary of 
region HI and IV for example. Although the turbulent spec- 
tra are described by similar spectral indices, the nonlinear 
processes are fundamentally different in region HI and IV. 
With this paper, we attempt to identify and understand the 
processes that lead to the spectra in HI and IV. It appears 
that the character of the nonlinear interaction is determined 
primarily by the orientation of turbulent velocity (V) and 
magnetic (B) field fluctuations with respect to each other, 
relative to the mean magnetic field. The spatially varying 
angular distribution of perpendicular velocity and magnetic 
field fiuctuations relative to the mean magnetic field was pre- 
dicted theoretically by Boldyrev (2006) and is a conjecture 
by Podesta et al (2008). It nevertheless remains to be seen 
how the orientation of turbulent V and B fluctuations govern 
the nonlinear spectral transfer of energy in solar wind tur- 
bulent plasma. In the context of MHD turbulence, Servidio 
et al (2008) show that the orientation between the velocity 
and magnetic fleld fluctuations plays a critical role in deplet- 
ing the nonlinear interactions that lead to the relaxation of 
MHD turbulence. 

Spectral anisotropy is another issue that is not yet prop- 
erly understood for characteristic modes with scales greater 
than kdi ~ 1 in Hall MHD plasma. In the case of MHD 
turbulence, the presence of a mean magnetic field leads to 
the asymmetric transfer of spectral energy alo ng and across 
the mean magnetic field l|Shebalin et al. 19831 ). An excellent 



analysis of anisotropic cascades in the framework of Hall 
MHD is presented by Ghosh & Goldstein (1997). They de- 
scribe how the degree of turbulent anisotropy in a Hall MHD 
plasma varies as the plasma beta (the ratio of magnetic to 
pressure energy) changes from greater than to less than 1 . A 
more general characterization of spectral anisotropy in the 
vicinity of modes kdi ~ 1 is still unclear. Specifically, the 
scale dependence of the turbulent anisotropy in the KAW 
regime and it's connection with that in the whistler regime 
is not yet established. Spectral anisotropy in the whistler 
regime was independently investigated by Shaikh & Zank 
(2003), Dastgeer et al (2000a, b). In this paper, we re- 
late these results to those associated with anisotropic KAW 
modes. 

To address the issues described above, we use time de- 
pendent, fully compressible three dimensional simulations of 
Hall MHD plasma in a triply periodic domain. This repre- 
sents a local or regional volume of the solar wind plasma. 
Note that the dynamics of length-scales associated with re- 
gion HI, i.e. corresponding to the KAW modes, cannot be 
described by the usual MHD models since we are inter- 
ested in characteristic frequencies smaller than an ion gyro 
frequency. At 1 AU, ion inertial len gth scales are smalle r 
than ion gyro radii in the solar wind l|Goldstein et al 19951 ). 
Plasma effects due to finite Larmor radii can readily be in- 
corporated in MHD models by introducing Hall terms to 
accommodate ion gyro scales up to scales as small as ion 
inertial length scales. In section 3, we describe the underly- 
ing Hall MHD model along with the interinsic assumptions 
and discuss the Hall MHD linear dispersion relation. Section 
4 addresses our 3D fluid simulations of the nonlinear Hall 
MHD equations. Our results suggest that the secondary iner- 
tial range spectrum has a form that is close to k"''^^ above 
the spectral break in the solar wind plasm, and mediated 
by the Hall terms in the short wavelength (in comparison 
with the ion skin-depth) KAW regime. Anisotropy in spec- 
tral cascade behavior is explored in section 5. We flnd that 
long length scale fluctuations in the kdi > 1 KAW regime ex- 
hibit a more anisotropic energy cascade compared to smaller 
scales. The dynamical alignment and angular distribution 
of turbulent velocity and magnetic fleld fluctuations is de- 
scribed in section 6. We find that characteristic turbulent 
fiucutations in the kdi > 1 regime relax towards a state of 
orthogonality such that the majority of turbulent scales con- 
tain fluctuations in which the velocity and magnetic flelds 
are nearly orthogonal , i.e. V _L B. Section 7 compares the 
KAW and whistler spectra. Finally, in section 8 we provide 
a summary and conclusions. 



3 HALL MHD SIMULATION MODEL 

Our 3D simulations are based on a two fluid nonlin- 
ear Hall MHD plasma model. The model assumes that 
the electrons are inertial- less, while the ions are inertial 
terishan fc Mahaian 2003 ). Hence, the electrons and ions 
have a differential drift, unlike the one fluid MHD model for 
which the electron and ion flow velocities are identical. Hall 
MHD description of magnetized plasma has previously been 
employed by a number of workers to investigate wave and 
turbulence processes in the context of solar wind plasma. 
In an excellent work, Sahraoui et al. (2007) extended the 
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ordinary MHD system to include spatial scales down to the 
ion skin depth or frequencies comparable to the ion gyrofre- 
quency in am incompressible limit. They further analyzed 
the differences in the incompressible Hall MHD and MHD 
models within the frame work of linear modes, their dis- 
persion and polarizations. Galtier (2006) developed a wave 
turbulence theory in the context of an incompressible Hall 
MHD system to examine the steepening of the magnetic 
fluctuation power law spectra in the solar wind plasma. Fur- 
thermore, Galtier and Buchlin (2007) have developed 3D 
dispersive Hall magnetohydrodynamics simulations within 
the paradigm of a highly turbulent shell model and demon- 
strated that the large-scale magnetic fluctuations are char- 
acterized by a fc^^''^-type spectrum that steepens at scales 
smaller than the ion inertial length di to k~''^^ . 

Here we start from the electron momentum equation, 

mefle + Ve • v) Ve = -VP^-eUe (e + ^Ve X (1) 

where rn^, n^, Vg, arc respectively mass, density, velocity 
and pressure of electrons, and E and B are the electric and 
magnetic fields. In the presence of low-frequency (compared 
with the electron gyrofrequency) electromagnetic fields, the 
electric force acting on the inertial-less electrons is balanced 
by the electron Lorentz force and the pressure gradient. This 
yields a more general form of Ohm's law than is typically 
used in the MHD description and is described as Hall MHD. 
We assume a quasineutral solar wind plasma where the den- 
sity of electrons (rte) and ions {rii) is nearly equal such that 
Ue ~ Ui = n. Tims in the inertial-less (Ohm's law) electron 
limit, the electron momentum equation yields the electric 
field as 



(2) 



E= -VPe - -Ve X B 

ne c 

The electric field arising from separation introduced by the 
inertial-less electron momentum equation can be substituted 
into the ion momentum equation 

mim + Vi • Vi = -VPi - em (b + i V, x B j , (3) 

where me, ne, Ve, Pc arc mass, density, velocity and ion pres- 
sure respectively, which then yields 

d 



1 



■ + Vi • Vi = -VP + ^ J X B. 



(4) 



Here we used p = rriin (solar wind plasma mass density), 
J = ne(Vi — Ve) (plasma current) and P = Pe + Pi (total 
plasma pressure). We next substitute the electric field into 
Faraday's law which yields 



^i- = -cV X E = V X (Ve X B). 



(5) 



On eliminating the electron fluid velocity from the current 
equation i.e. Ve = Vi — J/ne, we obtain 



m 

dt 



V X 



(V, 



X B 



J X B 



(6) 



In the inertial-less electron (me — > 0) limit the electron 
fluid does not influence the momentum of solar wind plasma 
directly except through the current. Since the electron fluid 
contributes to the electric field, plasma currents and the 
magnetic field are affected by electron oscillations. The com- 
bination of electron dynamics and ion motions distinguishes 



the Hall MHD model from its single fluid MHD counterpart. 
Thanks to the inclusion of electron dynamics, Hall MHD can 
describe solar wind plasma fluctuations that are associated 
with a finite ion Larmor radius and thus a characteristic 
plasma frequency is > u)ci- Because Hall MHD contains 
both ion and electron effects, there is a regime at which the 
dominance of one set of plasma fluctuations crosses over to 
be dominated by the other. This therefore introduces nat- 
urally an intrinsic scale length/timescale (frequency) that 
separates ion dominated behavior in the plasma from elec- 
tron dominated. It is the Hall term corresponding to the 
J X B term in Faraday's law that is primarily responsible 
for decoupling electron and ion motion on ion inertial length 
and ion cyclotron time scales (and introducing an intrinsic 
length scale). It is this feature that makes Hall MHD useful 
in describing dissipative solar wind processes when single 
fluid MHD is not applicable (the MHD model breaks down 
at a; > ujci)- Hall MHD allows us to study inertial range 
cascades beyond u) > uid, and can be extended to study 
dissipative heating processes where ion cyclotron waves are 
damped. The extreme limit of a fluid modeling applied to 
solar wind processes (even beyond the limit of the Hall MHD 
regime) is the use of an electron MHD model in which high 
frequency electron dynamics is treated by assuming station- 
ary ions that act to neutralize the plasma background. We 
discuss this regime at the end of this paper. 

The quasi- neutral plasma density (p), the velocity (V), 
the magnetic field (B) and the plasma pressure (P = Pe + 
Pi) can be described by the nonlinear Hall MHD equations. 
For the purpose of numerical simulations, it is convenient to 
express Hall MHD in the conservative form, 

f+V.Q = Q, (7) 
where, 



F = 



Q 



P 

pY 
B 



,Q 



pV 

PVV+^J+^P 



VB-BV 



-di(JB 
W 



BB 
BJ) 







fM(r,t)-hAtV'V-Fr?V(V- V) 
??V^B 



+ 



87r' 



and 

^=2^^ +(^ 
^4^^'^+(^^ + l^^><^- 

The suflix i is dropped from the ion fluid velocity Vi in the 
conservative form. 

The dynamical variables are functions of three space 
coordinates and time, i.e. {x,y,z,t) and arc normalized 
by typical length and time to = £o/Va in our simula- 
tions, Va = Bo/{A-kpq)^^^ the Alfven speed, such that V = 
£oV,d/dt = tod/dt,V = V/Va,B = B/VA(47rpo)'/', E = 
E/VA(47rpo)'/',P = P/poVlp = p/po,c = c/Va. I is the 
unit tensor. The parameters p, and r] represent the ion- 
electron viscous drag and magnetic field difFusivity, respec- 
tively. While the viscous drag modifies the dissipation in 
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the plasma momentum in a nonlinear fashion, the mag- 
netic diffusivity damps small scale magnetic field fluctu- 
ations linearly. The dimensionless parameter di(= di/lo, 
where di = c/ujpi is the ion skin depth and uipi is the ion 
plasma frequency) in Faraday's law identifies the Hall effect. 
The ion skin depth is a natural or an intrinsic length scale 
representative of the Hall MHD plasma model, and the Hall 
term plays an important role for high-frequency fluctuations 
with kdi ^ 1. It turns out that the Hall physics dominates 
the magnetoplasma dynamics when (l/p)J x B > V x B 
in the magnetic field induction equation. Finally, our non- 
linear Hall MHD model also includes the full energy equa- 
tion, rather than an adiabatic equation of state connecting 
the plasma pressure and the density. The use of the energy 
equation enables us to follow self-consistently the evolution 
of turbulent plasma heating resulting from the nonlinear cas- 
cading of energy. 

The rhs in the momentum equation corresponds to a 
forcing function (fA/(r, t)) that drives the plasma momentum 
at large length scales in our simulation. With the help of this 
function, we introduce energy in large scale eddies to sus- 
tain the magnetized turbulent interactions. In the absence 
of forcing, the turbulence decays freely. While the driving 
term modifies the plasma momentum, we conserve density 
(since we neglect photoionization and recombination). The 
large-scale random driving of turbulence can correspond to 
external forces or instabilities, for example fast and slow 
streams, merged interaction regions etc, in the solar wind, 
(or even supernova explosions and stellar winds in the ISM, 
etc). 

The linear Hall MHD equations admit the dispersion 
relation (|Brodin fc Stenflo lOQa l. 

2 

which has been written in a form to clearly identify the 
relationship of MHD to Hall MHD. The expression 

Dm{LO, k) = - u^k^iVl + V!) + kl^Vlv! 

is the MHD dispersion relation for fast and slow magne- 
tosonic modes. Here uj is the wave frequency, k(= kx -f kzi) 
the wave vector, the subscripts _L and z represent compo- 
nents across and along the external magnetic fields Bqz, and 
Va and Va are the Alfven and sound speeds respectively. 
Besides the fast and slow magnetosonic modes, KAWs, EM- 
ICA waves, and electron whistlers are identified by the re- 
maining three roots of the above dispersion relation. Warm 
Hall MHD plasma thus supports a great variety of waves 
of different wavelengths (comparable to the ion skin depth 
associated with the Hall drift, the ion sound gyroradius as- 
sociated with the electron pressure and the perpendicular 
ion inertia/ion polarization drift). 



4 NONLINEAR SIMULATION RESULTS 

To study the evolution of turbulent cascades in the solar 
wind turbulent plasma, we have developed a fully 3D com- 
pressible Hall MHD code. The three dimensional compu- 
tations are numerically expensive. But, with the advent of 
high speed vector and parallel distributed memory clusters, 
and efficient numerical algorithms such as those designed for 



Figure 2. Scaling of our three dimensional Hall MHD code with 
respect to number of processors. The spectral resolution is fixed 
at 60"^ in a cubic box of volume ir^ . By increasing the number of 
shared memory processors, the total computational time (mea- 
sured in CPU unit time; real time) to achieve steady state tur- 
bulent interactions can be reduced significantly. Our code clearly 
scales efficiently with an increasing number of processors. 



Message Passing Interface (MPI) libraries, it is now possi- 
ble to perform magnetofluid turbulence studies at substan- 
tially higher resolutions. Based on MPI libraries, we have de- 
veloped a three dimensional, time dependent, compressible, 
non-adiabatic, driven and fully parallelized Hall magnetohy- 
drodynamic (MHD) nonlinear code that runs efficiently on 
both distributed memory clusters like distributed-memory 
supercomputers or shared memory parallel computers. This 
allows us to achieve a very high resolution in Fourier spec- 
tral space. The code is scalable and transportable to differ- 
ent cluster machines. The spatial descretization employs a 
conservative and parallelized descretization of Fourier har- 
monics, and the temporal integration uses a Runge Kutta 
(RK) 4th order method. The boundary conditions are peri- 
odic along the x, y and z directions in the local rectangular 
region of the solar wind plasma. The MHD counterpart of 
this code was used in several other studies by us, see e.g . 
(|Shaikh fc Zank 20061 : IShaikh fc Zank 20071 : IShaikh 2009"al ) . 

Our code treats the solar wind plasma fluctuations as 
statistically isotropic, locally anisotropic, homogeneous and 
random. Such a representation is further consistent with 
ACE spacecraft measurements (Smith et al 2006). The nu- 
merical algorithm accurately describes the physical variables 
(density, temperature, magnetic field, velocity field, pressure 
etc) in our code and they are also less dissipative. Because 
of the latter, nonlinear mode coupling interactions preserve 
ideal rugged invariants of fiuid fiows, unlike finite difference 
or finite volume methods. The conservation of ideal invari- 
ants (energy, enstrophy, magnetic potential, helicity) in in- 
ertial range turbulence is an extremely important feature 
because these quantities describe the cascade of energy in 
the inertial regime, where turbulence is, in principle, free 
from large-scale forcing as well as small scale dissipation. 
Damping of plasma fiuctuations may nonetheless occur as a 
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result of intrinsic non-ideal effects such those introduced by 
the finite Larmor radius. 

To test the scalability of our code, we have performed 
simulations by increasing the number of processors for a 
fixed number of modes. Our scaling results are depicted in 
Fig. Clearly, by increasing the number of shared memory 
processors, the computational time to achieve steady state 
turbulent interactions can be reduced significantly. 

The turbulent fluctuations are initialized by using a 
uniform isotropic random spectral distribution of Fourier 
modes concentrated in a smaller band of lower wavenum- 
bers (fc < 0.1 kmax)- While spectral amplitudes of the fluc- 
tuations are random for each Fourier coefficient, it follows 
a certain initial spectral distribution proportional to k~°', 
where a is the initial spectral index. The spectral distri- 
bution set up in this manner initializes random scale tur- 
bulent fluctuations. We note that a constant magnetic field 
is included along the z direction (i.e. Bo = Boz) to ac- 
commodate the large scale (or the background solar wind) 
magnetic field. The plasma beta, ratio of plasma pressure 
to magnetic field energy, is close to unity in our simula- 
tion /3 ~ 1, as typically observed in solar wind (Goldstein 
et al 1995; Smith et al 2006). Turbulent fiuctuations in our 
3D Hall MHD simulations are driven either at the lowest 
Fourier modes or evolve freely under the infiuence of the 
self-consistent dynamics described by the set of Equation 
([7|. The inertial range spectral cascades in either cases lead 
to nearly identical inertial-range turbulent spectra. We have 
further carried out 3D simulations for a range of various pa- 
rameters and spectral distributions to ensure the validity, as 
well as the consistency, of our codes and the physical results. 
The simulation parameters are; spectral resolution is 128^, 
77 = /i = I0''\f3 = 1.0, kdi ~ 0.1 - 10, Lx = Ly = L:, = 2tt. 
The nonlinear coupling of velocity and magnetic field fluctu- 
ations, amidst density perturbations, excites high-frequency 
and short wavelength (by the ui/uJci effect) compressional 
dispersive KAWs. While the nonlinear mode coupling in- 
teractions must influence local turbulent fluctuations in the 
inertial range, their role in the spectral energy cascade is not 
yet fully understood. We come back to this issue below. 

The nonlinear spectral cascade in the modified KAW 
regime leads to a secondary inertial range in the vicinity of 
kdi — 1, where the turbulent magnetic and velocity fluctu- 
ations form spectra close to k~''^^. This is displayed in Fig. 
((21, which also shows that for length scales larger than the 
ion thermal gyroradius lead, an MHD inertial range spec- 
trum close to k~^^^ is formed. Figure (|3} results from aver- 
aging the spectra over ten simulation runs that are initial- 
ized with different random numbers. The spectra obtained 
in each simulation further consist of 20 different spectra that 
are averaged over the simulation period (15 — 20) Iq/Va- The 
resultant spectrum is further processed using a higher order 
polynomial flt leading thereby to a well-behaved omnidi- 
rectional power law. The characteristic turbulent spectrum 
in the KAW regime is steeper than that of the MHD in- 
ertial range. Identifying the onset of the secondary inertial 
range has been the subject of some debate because of the 
presence of multiple processes in the KAW regime that can 
mediate the spectral transfer of energy. These processes in- 
clude, for instance, the dispersion and damping of EMICA 
waves, turbulent dissipation, etc. In the context of our 3D 
Hall MHD simulations, the observed A;"^''^ scaling in the 



Landscape/two column figure to go here 

Figure 3. Inertial range turbulent spectra for magnetic and ve- 
locity field fluctuations. The fluctuations closely follow respec- 
tively fc"^/^ and k'"^/^ scalings in the kdi < 1 and kdi > 1 KAW 
regimes, di = 0.05 and 1.0 respectively in the kdi < 1 and kdi > 1 
regimes. The dash-dot straight lines correspond to a k~^^'^ and a 
fc— ■^/S power law. 

turbulent plasma can be understood from effect of the Hall 
term on the energy cascade process. The time-scale associ- 
ated with Hall MHD is th, which is smaller than the MHD 
characteristic time scale tmhd. Nonlinear cascades in MHD 
turbulence are typically governed by 

1 

Tmhd ~ -; — , 
kvk 

where Vk — |V(k,t)| is the velocity field in fc-space. By con- 
trast, spectral transfer of turbulent energy in the Hall MHD 
plasma has a typical timescale of 

1 



The energy transfer rate in the KAW regime is, therefore, 

(8) 

On assuming a t urb ulent equipartition relation 
(|Alexandrova et al 20071 : lAlexandrova et al 20081 ) be- 
tween the velocity and magnetic fields ~ v1, and using 
Th, the energy transfer rate, Eq. (|8]) reduces to 

e ~ k'^Bt 

The inertial range modes in our 3D Hall MHD 
simulations exhibit near equipa rtition, a result that 
is consistent with ot her work (|Alexandrova et "aT 20071 : 
lAlexandrova et al 20081 ) . This is shown in Fig. Q. It is note- 
worthy that equipartition was assumed by Alexandrova et 
al (2007) and Alexandrova et al (2008) and is not directly 
observed in the solar wind plasma in the kdi > 1 regime. 
Our simulation results shown in Fig. Q thus provide theo- 
retical support to arguments that i nvoke turbulent equipar- 
tition in the MHD/KAW range the (|Alexandrova et al 20071 : 
lAlexandrova et al 20081 ). Fig. Q shows turbulent equipar- 
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tition where spectrally averaged (i.e. averaged over entire k- 
spectrum) fluctuations exhibit turbulent equipartition dur- 
ing the nonlinear evolution. In Fig. ([4]), the ratio of mag- 
netic and velocity field energy remains close to unity i.e. 
|-BfcP/|i'fcP ~ 1. Fig. ^ describes spectral distribution of 
\5vf/\SB\^ as a function of k. This fi gure also shows that 
the ratio |(5v|^/|(5-B|^ stays more or less close to the unity for 
the inertial range turbulent scales {kdi ^ 1) in our simula- 
tions. 

In the kdi < 1 Alfvenic regime, turbulent equiparti- 
tion i n the solar wind plasma remains a subject of de- 
bate JPodesta et al 200"! : ITu fc Marsch 19931 : lYokoi 20061 : 
lYokoi fc Hamba 2007 ). Turbulent equipartition observed in 
the solar wind plasma in the kdi < 1 Alfvenic regime 
is not exact and vl/B^ ^ 1. A precise mechanism lead- 
ing to the non-exact turbulent equipartition is unclear in 
the solar wind plasma (in the kdi < 1 Alfvenic regime). 
The deviation from equipartitio n may resu lt from inho- 
mogeneity of the magnetic field (|Yokoi 20061 ). possible dif- 
feren ces in the dissipation of kinetic and magnetic ener- 
gies (|Yokoi fc Ham ba 2007). or compressibility of the den- 
sity fluctuations (Podcsta ct al 2008). Our objective how- 
ever is not to address the mechanism that lead to un- 
equal turbulent equipartition in the solar wind plasma 
in the kdi < 1 Alfvenic regime. The readers can re- 
fer to the work of /Podesta et al 20081 : ITu fc Marsch 19931 : 
lYokoi 20061 : lYokoi fc Hamba 2007I T for a more detailed anal- 
ysis of turbulent equipartition in the kdi < 1 Alfvenic regime 
in the solar wind plasma. We simply use the equipartition 
relationship here to derive Kolmogorov-like energy spectrum 
in Hall MHD plasma. 

To derive the expression for the energy spectrum 
describing the energy cascade in H all MHD plasma , 
we apply Kolmogorov's phenome nology (|Kolmogorov 194ll : 
iKraichnan 19651 : llroshnikov 19631 ) in which the energy cas- 
cade in the inertial range is local and d epends on Fourier 
modes and the e nergy dissipation rates ([Kolmogorov 194ll : 
IKraichnan 19651 : llroshnikov 1963t ). We obtain 

k-'Bl 



{Bierk" 



On equating indices, we obtain a = 2/3 and /3 = 
results in an energy spectrum of the form 



-7/3. This 



which is consistent with our 3D Hall MHD simulations [see 
Fig. (|3J] and solar wind observations (Goldstein et al 1995; 
Leamon et al 1999, Bale et al. 2005, Sahraoui et al 2009, 
Alexandrova et al 2007, 2008). On the other hand, the use 
of tmhd in estimating the energy dissipation rates recovers 
the k~^^^ inertial range MHD spectrum. This suggests that 
the Hall effect may be responsible for the spectral steepening 
in the solar wind plasma fluctuations in the kdi > 1 regime. 

It is worth mentioning that Fig. Q represents spectral 
average of the entire mode spectrum. On the other hand, 
spectra of \5B{k)\^ /\5V{k)\^ and \B{k)\^ /\V{k)\^ in steady 
state show that turbulent equipartition in Hall MHD devi- 
ates for higher fc's (Shaikh fc Shukla 2009). Interestingly, 
nonlinear fluid simulations in whistler wave regime (Dast- 
geer et al 2000a, Shaikh 2009) corresponding to region IV 
in Fig (1) show that magnetic and velocity fields in whistler 
modes tend to establish turbulent equipartition. This is es- 
sentially a regime where the characteristic length scales are 



Figure 4. Volume averaged spectra of velocity and magnetic field 
fluctuations exhibit turbulent equipartition as a function of time. 



Figure 5. Spectral distribution of |(5i)p/|<5Bp as a function of 
k. Figure shows that the ratio |(5i)p/|<5Bp stays more or less close 
to the unity. 



bigger than electron inertial length, i.e. kde < 1. In this 
regime, dispersive wave effects dominate the cascade pro- 
cesses by establishing turbulent equipartion amongst modes. 
In the fcde > 1 regime (region V in Fig (1)), Dastgeer et al 
(2000b) find that turbulent equipartion does not hold for 
short scales in EMHD turbulence. The latter behaves like 
hydrdynamic regime of fluid turbulence where total energy 
is dominated entirely by kinetic energy, whereas magnetic 
energy becomes sub-dominant. 



5 ANISOTROPIC CASCADES 

Here we discuss the degree of anisotropy in the turbulence, 
introduced by the presence of a large scale magnetic field in 
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3D Hall MHD. A mean magnetic field in the z-direction is 
included in our simulation, this meant to mimic the large 
scale solar wind background magnetic field. The small scale 
fiuctuations are infiuenced by the presence of the large scale 
mean magnetic field and tend to evolve in an anisotropic 
manner in the sense that the turbulent cascade along and 
across the mean magnetic field behave difi'erently. In 3D 
turbulence, fluctuations in a plane orthogonal to the mean 
magnetic field remains isotropic, whereas those aligned in 
the direction of the external magnetic field are affected. The 
latter is establised for /? < 1 and /3 ~ 1 by Zank & Matthaeus 
(1990, 1993) in the context of a nearly incompressible MHD 
theory. Thus, turbulent anisotropy in the inertial range spec- 
trum corresponds to the preferential transfer of spectral en- 
ergy that feeds perpendicular modes k± while the parallel 
(fc||) cascades are suppressed. The anisotropy in the initially 
isotropic turbulent spectrum is triggered essentially by back- 
ground anisotropic gradients that nonlinearly migrate spec- 
tral energy in a particular direction. To measure the degree 
of anisotropic cascades, we employ the following diagnostics 
to monitor the evolution of the k± mode in time. The av- 
eraged k± mode is determined by averaging over the entire 
turbulent spectrum weighted by thus 



(fcx(t)> 



Here (• ■ •) represents an average over the entire Fourier spec- 
trum, kj_ = yfci + fcf and Q represents any oi B,V, p, V x B 
and V X V. Similarly, the evolution of the fe|| mode is deter- 
mined by the following relation. 



(fe|iW> 



TJfciiQ(fc>^)P 



It is clear from these expressions that the {k±{t)) and (fcy (t)) 
modes exhibit isotropy when {k±{t)) ~ Any devia- 

tion from this equality corresponds to spectral anisotropy. 
We follow the evolution of the {k±{t)) and modes 
in our simulations. Our simulation results describing the 
evolution of {k±) and (fc||(t)) modes are shown in Fig. 
((6|. It is evident from Fig. (|6]) that the initially isotropic 
modes {k±{t)) ~ gradually evolve towards a highly 

anisotropic state in that spectral transfer preferentially oc- 
curs in the {k±{t)) mode, and is suppressed in mode. 
Consequently, spectral transfer in the {k±{t)) mode domi- 
nates the nonlinear evolution of fluctuations in Hall MHD, 
and mode structures become elongated along the mean mag- 
netic fleld or 2-direction. Hence nonlinear interactions led by 
the nonlinear terms in the presence of background gradients 
lead to anisotropic turbulent cascades in the inertial range 
turbulent spectra. 

While the spectral transfer of energy differs in the k± 
and fcy modes, the 3D volume averaged turbulent spectrum 
follows a k~'^^^ power law, as shown in Fig. ((Sj. The steep- 
ness of the observed spectrum can be ascribed to the co- 
existence of partially anisotropic flows and turbulent fiuc- 
tuations in steady state Hall MHD turbulence. Fig. (|6]) il- 
lustrates anisotropy corresponding to an averaged k mode 
but the anisotropy exhibited by the small and large scale 
B and v fluctuations is not distinctively clear, nor is the 
degree of anisotropy in B and v flelds clear from Fig. (|6]). 



Figure 6. Evolution of {k_i_{t)) and (fc||(t)) as a measure of 
anisotropy in Hall MHD. Initially, {k^{t = 0)> = k\\(t = 0). 
As time evolves, anisotropy in spectral transfer progressively de- 
velops such that {k±{t)) > {fc||(t)). Hence the presence of a large 
scale magnetic field suppresses the turbulent cascade along the 
direction of the field, while the perpendicular cascade remains 
nearly unafi^ected. 



Figure 7. Spectrum of kj_ and fcy as a function of k. The large 
scale inertial range turbulent fluctuations are more anisotropic as 
compared to the smaller ones. 



The scale dependece of turbulent anisotropy is described in 
Fig. ((TJ. Fig. (O shows discrepancy in k± and fcy is promi- 
nent at the smaller k's. This essentially means that the large 
scale turbulent fluctuations are more anisotropic than the 
smaller ones in a regime where characteristic length scales 
are smaller than di i.e. kdi > 1. It further appears from 
Fig. (O that the smaller scales in the kdi > 1 are virtually 
unaffected by anisotropic kinetic Alfven waves that propa- 
gate along the externally imposed mean magnetic fleld Bq. 
Turbulent fluctuations with small characteristic scales in the 
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Figure 8. Evolution of the angular anisotropic distribution in 
Hall MHD turbulent plasma fluctuations. Shown in the figure 
is the angular evolution of anisotropy associated with the rota- 
tional velocity and magnetic field fiuctuations. The initial angu- 
lar distribution is identical. As the simulation progresses, small 
scale velocity field fluctuations tend to become more anisotropic 
than their magnetic field counterpart. The implication is that the 
smaller scale magnetic field fluctuations are less anisotropic than 
the velocity fleld fluctuations. 



kdi > 1 regime of Hall MHD are not affected by the mean 
magnetic field or kinetic Alfven waves. This leads us to con- 
jecture that small scale turbulence in the kdi > 1 regime 
behaves essentially hydrodynamically i.e. as eddies indepen- 
dent of the mean magnetic field or coUisionless magnetized 
waves. Thus large and smaller turbulent length scales evolve 
differently in the kdi > 1 regime of Hall MHD turbulence. 

Consider now the anisotropic evolution of B and V tur- 
bulent fiuctuations in the Hall MHD regime of solar wind 
plasma, i.e. the kdi > 1 wavenumber band. Specifically we 
investigate the degree of anisotropy associated with the non- 
linear transfer of turbulent energy in these fiuctuations. The 
results of our 3D fully compressible Hall MHD plasma sim- 
ulations are plotted in Fig. ((S]). In Fig. we quantify 
the turbulent anisotropy associated with small scale turbu- 
lent fluctuations that correspond to the rotational B and V 
fields to investigate which is more anisotropic in the pres- 
ence of the mean magnetic field. The evolution of turbulent 
anisotropy, corresponding to small scale rotational B and V 
fiuctuations, is determined from the following relations, 



■ tan 



tan 



^JfcxpKkxB(k,f)p 

^Jfc||P|ikxB(k,t)|2 

EJfc^riikx V(k,t)p 

^Jfc||P|jkX V(k,i)|2 



1/2 



1/2 



(9) 



(10) 



For isotropy, k± ~ fcy , hence the angle of ansiotropy is 
close to S ~ 45° . Any deviation from the isotropic spectrum 
can be seen if 6* / 45°. It is also clear from Eqs. ^ k. (fTO)) 
that 9 > 45° corresponds to the inertial range anisotropic 
turbulent spectrum that is dominated by the cascades in the 
perpendicular direction such that k±_ > fcii whereas 9 < 45° 



Figure 9. Evolution of the degree of alignment in Hall MHD fluc- 
tuations. A progressive decrease in the angle of alignment from 
90° indicates the eventual weakening of V X B nonlinear interac- 
tions in the Alfvenic cascade regime of solar wind turbulence. 



describes k± < fcy cascades. The evolution of the angular 
anisotropy associated with small scale inertial range turbu- 
lent fiuctuations described by Eqs. (O & (|lUp is depicted 
in Fig. ((8)1. Several important points emerge from Fig. (|8]). 
Firstly, spectral anisotropy associated with small scale in- 
ertial range turbulent fiuctuations in the velocity field is 
more pronounced than that in the magnetic field fiuctua- 
tions at large times. The angular evolution associated with 
the anisotropic velocity field deviates markedly from the 
isotropic angle — 45°. On the other hand, magnetic field 
fiuctuations show a smaller degree of anisotropy compared 
to the velocity field fiuctuations. It then appears, in view of 
Fig. ([7)l, that the velocity field spectrum is dominated by 
relatively large scale turbulent fiuctuations as compared to 
the magnetic field. 



6 DYNAMICAL ALIGNMENT OF 
TURBULENT FLUCTUATIONS 

The nature and strength of nonlinear interactions resides 
critically with their orientation with respect to each other 
in the presence of an externally imposed or self-consistently 
generated large scale magnetic field. For instance, the domi- 
nant nonlinear interactions in a Hall MHD fiuid are governed 
by the Lorentz force that contains a nonlinear term corre- 
sponding to the V X B nonlinearity. Hence the orientation of 
the velocity field relative to the magnetic field fiuctuations 
plays a critical role in determining the strength of nonlin- 
ear interactions. This clearly means that nonlinear interac- 
tions are dominated by those fiuctuations that possess the 
velocity field aligned perfectly orthogonal (at a 90° angle) 
to the magnetic field, i.e. V ± B. The obliqueness, intro- 
duced primarily by deviation from the orthogonality, tends 
to weaken the strength of nonlinear interactions mediated 
by the V x B nonlinearity. The V x B nonlinearity disap- 
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Figure 10. Evolution of the degree of alignment in the kdi > 1 
regime of Hall MHD. Small scale fluctuations possess orthogonal 
velocity and magnetic fields and hence the angle of alignment 
stays close to 90°. 



characteristic turbulent length scales in the kdi > 1 regime. 
The angle of alignment for the smaller scales correspond- 
ing to the kdi > 1 regime is shown in Fig. (I10|l . Significant 
differences are apparent in the angle of alignments associ- 
ated with the large and small scales Figs (|9]) & (fTO)) . It ap- 
pears from our simulations that the small scale fluctuations 
{kdi > 1) are nearly orthogonal as seen in Fig (10). By con- 
trast, the large scale fluctuations {kdi < 1) in Fig (9) show 
a signiflcant departure from the orthogonality. 

It is interesting to note that characteristic length scales 
that are large compared to the ion inertial skin depth 
{kdi < 1) tend to deviate from orthogonality (a tendency to- 
wards dynamic alignment). Furthermore, these are the scales 
at which anisotropic cascades are dominated by the higher 
rate of spectral transfer of energy in the k± modes. This sug- 
gests the possibility that anisotropic cascades are related to 
the orientation of the velocity field relative to the magnetic 
field fiuctuations. Our simulation results described in Figs 
(0, (O, (O 1^ (|10|> hint that the predominance of spectral 
anisotropy at large scales and the obliqueness are related. 
A heuristic and theoretical analysis of this relationship is 
beyond the scope of this paper. 



pears (i.e. V |[ B = 0) for fluctuations that carry parallel or 
anti-parallel velocity and magnetic fleld fluctuations. 

To understand the strength of the nonlinear interactions 
in Hall MHD solar wind plasma, we determine the degree of 
alignment of the velocity and magnetic field fluctuations by 
defining the following alignment parameter (Podesta et al 
2008) that spans the entire fc-spectrum in both the kdi > 1 
(HaU MHD) and kdi < 1 (usual MHD) regimes. 



e(t) = cos" 



Ei.Vk(t)-Bi.(t) 
Ek|Vk(t)||Bk(f)| 



The summation is determined from the modes by summing 
over the entire spectrum. In this sense, the alignment pa- 
rameter depicts an average evolution of the alignment of ve- 
locity relative to the magnetic fleld fluctuations. Note care- 
fully that this alignment can vary locally from smaller to 
larger scales, but the averaging (i.e. summing over the entire 
spectrum) rules out any such possibility in our simulations. 
Nonetheless, Q deflned as above enables us to quantitatively 
measure the average alignment of the magnetic and velocity 
field fluctuations while the nonlinear interactions evolve in 
a turbulent solar wind plasma. 

Unlike Podesta et al. (2008), where a probability distri- 
bution of angular orientation corresponding to each mode is 
determined, we follow the total (i.e. spectrally volume av- 
eraged) evolution of the angle of alignment for the velocity 
and magnetic fleld fluctuations in both the Alfvenic cascade 
(i.e. MHD) and kinetic Alfven wave (i.e. Hall MHD) regimes. 
Our results are plotted respectively in Figs © & pop . The 
Alfvenic cascade regime of MHD turbulence kdi < 1 in the 
solar wind plasma possesses relatively large scales {kdi < 1) 
in which the velocity and magnetic fleld fluctuations are 
observed to be somewhat obliquely aligned. Hence our sim- 
ulations show that the angle of alignment evolves towards 
< 90°, as depicted in Fig. ([9}. Hence the strength of the 
nonlinear interactions corresponding to the V x B nonlin- 
earity is relatively weak. This result is to be contrasted with 



7 COMPARISON BETWEEN KAW AND 
WHISTLER SPECTRA 

It is interesting to note the similarities and differ- 
ences between the characteristic Hall MHD and whistler 
wave turbulence spectra. Here we discuss properties 
of whistler turbulence from the previ o us works by 



or wnistier turpuience irom tne previ o us worKs py 
dShaikh fc Zank 20031: [shaikh fc Zank 20051 : IShaikh 2009"al : 
IShaikh 2009bl : iBiskamp et al 19961 )" and compare them with 
our present simulation results. 

While the two spectra correspond to essentially different 
frequencies as depicted schematically in Fig. ((1} , they are de- 
scribed by nearly identical spectral power laws in the vicinity 
of high frequency KAW modes (see the interface of boundary 
between HI and IV in Fig. ([T}). Thus, the small scale KAW 
in the kdi > 1 regime and relatively large scale whistler fluc- 
tuations exhibit a k^''^^ spectrum. Although the slope of the 
inertial range spectrum for these two modes is identical, the 
nonlinear physical processes yielding the spectra are distinc- 
tively different. Whistler modes are essentially governed by 
the motion of electrons in a static neutralizing ion back- 
ground. The characteristic frequencies therefore reside be- 
tween the ion and electron gyro frequencies {uod ^ w <C Wce) 
and characteristic length scales lie between the ion (di) and 
electron (de) inertial lengths i.e., di = c/u-pi < £ < de = 
c/ujpe, where cOpi and uipe are respectively the ion and elec- 
tron plasma frequency. By virtue of d^, there exist two in- 
ertial ranges that correspond to smaller kde > 1 and larger 
fcde < 1 lengthscales. Correspondingly, forward cascade tur- 
bulent spectra in these regimes exhibit k~^^^ and k~''^^ re- 
spectively as shown in schematic Fig. The whistler spec- 
trum k~'^^^ is produced essentially by fluctuations in elec- 
tron fluid while ions are at rest. This is de monstrated in 2D 
( Biskamp et al 19961 : IShaikh fc Zank 20051 ) as well as in 30 
(IShaikh 2009bl ) simulations. As described in (|Shaikh 2009bl ). 
the whistler waves cascade the inertial range turbulent en- 
ergy typically through the convective electron fluid velocity 
We ~ V X B. Thus the typical velocity of the magnetic fleld 
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eddy Be with a scale size £ can be represented by Ve — Be/£. 
The eddy turn-over time is then given by 



'^whit 



(11) 



It is this time scale that characterizes the nonlinear spectral 
transfer of energy in fully developed whistler wave turbu- 
lence. Kolmogorov-like dimensional arguments further yield 
a k~'^^^ spectrum in the kde < 1 characteristic turbulent 
regime, as described in Fig. ([T}. While the inertial range 
nonlinear cascades are determined essentially by the eddy 
turn over or spectral transfer time scale, the characteristic 
length scales longer and less than de are significantly influ- 
enced by the whistler in teraction time scales in a disparate 
manner (|Shaikh 2009bh . By contrast, as described above, 
the KAW spectrum mediated by Hall MHD processes in the 
kdi > 1 regime is produced by the combined motion of elec- 
tron and ion fluids amidst density fluctuations that evolve 
on dispersive KAW time and length scales. The KAW fluctu- 
ations relax towards a, k~''^^ spectrum in the kdi > 1 regime 
essentially by means of nonlinear mode coupling interactions 
in which both the electron and ion fluid perturbations par- 
ticipate on an equal footing. The most notable difference 
in the energy cascade processes associated with the KAW 
and whistler modes thus emerges from Equations ((8)) & (fTTj) 
where the energy transfer rates are determined typically by 
different nonlinear fluid velocities. 



8 SUMMARY 

In summary, we have investigated the nonlinear and tur- 
bulent behavior of a two fluid, compressible, three dimen- 
sional Hall MHD model. Hall MHD model is relevant to 
describe solar wind plasma and magnetospheric and labora- 
tory plasma environments. As exhibited by the solar wind 
plasma, multiple scale characteristic fluctuations lead to a 
complex power spectral density (or power spectrum) span- 
ning very low frequencies, corresponding to a flatter 
spectrum, followed by the usual "Alfvenic" spectrum cor- 
responding to a k~^^^ spectrum. How the Alfvenic cas- 
cade is terminated in the solar wind plasma, thus intro- 
ducing a spectral discontinuity near the ion gyro frequency, 
continues to challenge our understanding of nonlinear so- 
lar wind turbulent processes near the kdi ~ 1 band of 
wavenumbers. The slope of inertial range beyond the kdi ~ 1 
wavenumbers varies between -2 and -5 (Smith et al 1990, 
Goldstein et al 1994, Leamon et al 1999), depending upon 
what determines the underlying nonlinear interactions. Sev- 
eral mechanisms have been proposed to describe the spec- 
tral discontinuity of turbulent energy transfer in the solar 
wind plasma; these include kinetic Alfven waves (KAWs) 
(Hasegawa 1976), electromagneti c ion-cyclotron- Alfven 
(EMICA) waves |Garv et al 20081 : IWu fc Yoon 20071 ). non- 
linear HaU MHD interactions (Alexandrova et al 2007, 2008; 
Shaikh & Shukla 2008, 2008a), and suppression of Alfven 
fluctuations by proton cyclotron damping at intermediate 
wavenumbers (Stawicki et al 2001), so exciting whistler 
waves which are dispersive unlike Alfven waves. More com- 
plex interactions may be expected due to the presence 
of highly obliquely propagating un-d amped kinetic Alf ven 
wave (with to <^ ujd) as proposed by (|Howes et al. 20081 ). 



In this paper, our simulation results make interesting 
non trivial connections between regions that correspond to 
distinct wavenumbers in the composite spectra sketched in 
Fig. We find that the spectra exhibit a break, consis- 
tent with the spacecraft observations (Smith et al 1990, 
Goldstein et al 1994, Leamon et al 1999), in the vicinity 
of the kdi ~ 1 wavenumbers. The characteristic turbu- 
lent modes below these wavenumbers exhibit a k~^''^ like 
MHD spectrum. By contrast, turbulent fluctuations scale as 
k~'^^^ in the inertial range characterized by kdi > 1- The 
spectral index —7/3 is mediated essentially by equiparti- 
tion processes between the turbulent velocity and magnetic 
field fluctuations. Note however that turbulent equiparti- 
tion, for instance \Bkf — [VfcP, is not directly observed in 
the spacecraft data. Small scale inertial range fluctuations 
in the kdi > 1 regime tend to establish turbulent equipar- 
tition amongst the characteristic modes in our simulations 
The inertial range k~''^^ spectrum is a unique feature of 
nonlinearly interacting multi-scale electromagnetic fluctua- 
tions that follows from the turbulent equipartition of high 
frequency kinetic Alfven modes. We further find a consider- 
able departure from turbulent equipartition for the higher k 
modes. The departure from equipartition in our simulations 
can be ascribed essentially to compressibility, which we find 
to be enhanced for the higher k modes in the inertial range 
KAW spectra. 

In the presence of a large scale mean background 
magnetic field, small scale turbulent fluctuations exhibit 
anisotropic cascades. We find that the long length scales 
in the kdi > 1 KAW regime are more anisotropic compared 
to the shorter scales. Dynamical alignment and angular dis- 
tribution of turbulence velocity and magnetic field fluctu- 
ations is found to play a critical role in determining the 
degree of nonlinear interactions. We find that characteristic 
turbulent flucutations in the kdi > 1 regime relax towards 
orthogonality, so that most of turbulent scale fluctuations 
have velocity and magnetic fields that are nearly orthogo- 
nal , i.e. V _L B. For large scale fluctuations corresponding 
to the MHD regime, magnetic and velocity flelds are not 
perfectly orthogonal, being instead on average nearly 70° 
to each other. By contrast, small scale fluctuations in the 
kdi > 1 KAW regime exhibit nearly perfect orthogonality 
in that the average magnetic and velocity fields make an 
angle of nearly 90° with respect to each other. We finally 
noted the similarities and differences between the KAW and 
whistler spectra. Interestingly, the two modes exhibit a sim- 
ilar fc"^/'' spectrum in region III and IV of Fig. ((TJ that cor- 
responds respectively to the kdi > 1 and kde < 1 regimes 
in wavenumber space. While the KAW spectrum is domi- 
nated by the combined motion of ions and electrons, the 
whistler wave spectrum is governed predominantly by elec- 
tron motion only in the presence of a static neutralizing ion 
background. 

The spectral properties of nonlinear Hall MHD are par- 
ticularly relevant for understanding the observed solar wind 
and heliospheric turbulence. Hall MHD may also be use- 
ful for understanding multi-scale electromagnetic fluctua- 
tions and m agnetic field reconnect i on in the Earth's mag- 
netosphere (|Matthaeus et al 2005|; Phan et al 20061) and 
in laboratory plasmas ( Egcdal c t al 20071 : Qno et al 19961 : 
IHsu et al 2"oO0l : ICarter"et al 20061 ). 
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